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Symmetric Z2 spin liquids and their neighboring phases on triangular lattice 


Yuan-Ming Lu^ 

^Department of Physics, The Ohio State University, Columbus, OH 43210, USA 

Motivated by recent numerical discovery of a gapped spin liquid phase in spin-1/2 triangular- 
lattice J 1 -J 2 Heisenberg model, we classify symmetric Z 2 spin liquids on triangular lattice in the 
Abrikosov-fermion representation. We find 20 phases with distinct spinon symmetry quantum num¬ 
bers, 8 of which have their counterparts in the Schwinger-boson representation. Among them we 
identify 2 promising candidates and #20), which can realize a gapped Z 2 spin liquid with up 
to next nearest neighbor mean-field amplitudes. We analyze their neighboring magnetic orders and 
valence bond solid patterns, and hnd one state (#20) that is connected to 120-degree Neel order 
by a continuous quantum phase transition. We also identify gapped nematic Z 2 spin liquids in the 
neighborhood of the symmetric states and find 3 promising candidates (#1, #6 and #20). 

PACS numbers: 71.27.+a, 75.10.Kt 


I. INTRODUCTION 

Thanks to advances in numerical methods and in¬ 
creasing computational power, recently a lot of progress 
has been made in the pursuit of quantum spin 
liquidj^ (QSLs) in various frustrated quantum spin 
models^H^. In contrast to classical paramagnets driven 
by thermal fluctuations, QSLs are symmetric zero- 
temperature phases featured by c ollective excitations 
which obey fractio nal st atisticjl^lfl^and long-range quan¬ 
tum ent anglemenlP^li^ therein. Most recently using the 
density matrix renor maliz ation group (DMRG) method, 
two numerical studie^^*!^ of the spin—1/2 Heisenberg Ji- 
J 2 model on triangular lattice independentl y rep orted 
strong evidence for a gapped spin liquid phase^E^ in the 
finite range of J 2 IJ 1 , sandwiched by a noncollinear 120- 
degree Neel order at smaller J 2 /J 1 and a collinear stripy 
Neel order at larger J 2 I Ji- This is also supported by ear¬ 
lier variational Monte CarlcPti^ and othePHm studies 
on the same model. Due to its similarity to the gapped 
Z 2 s pin li quid phase discovered in kagome Heisenberg 
modeP*^, this discovery motivates us to look for sym¬ 
metric Z 2 spin liquid phases on triangular lattice as pos¬ 
sible candidates for the one realized in J 1 -J 2 model. 

A symmetric Z 2 spin liquid hosts 3 types of fraction¬ 
alized excitations: bosonic spinon b, bosonic vison v and 
fermionic spinon /. Each of them obeys mutual semionic 
statistic^Sl with the other two. Unlike usual ordered 
phases characterized by their broken symmetries, sym¬ 
metric spin liquids are characterized by the symmetry 
quantum numbers of their fractional excitations. In par¬ 
ticular, each quasiparticle can carry a fraction of the 
unit symmetry quantum number carried by a funda¬ 
mental particle (such as a magnon in frustrated mag¬ 
nets), these phe nomen a are generally coined “symmetry 

fractionalization’*2^f24|^ Take Z 2 spin liquids for instance, 

each spinon carries spin-1/2, half of the spin quantum 
number carried by each magnon. Taken into account of 
time reversal, spin rotation and space group symmetries. 


there are many different symmetric Z 2 spin liquids with 
distinct qusiparticle symmetry quantum numbers. In¬ 
stead of attempting to fully classifjJ^MUl possible Z 2 
spin liquids on triangular lattice, in this work we restrict 
ourselves to the Abrikosov-fermion representatio]j22HSl| 
which is a convenient way to construct variational wave- 
functions of Z 2 spin liquids. Among 20 distinct symmet¬ 
ric states in Abrikosov-fermion representation, we iden¬ 
tify two promising candidates (#1 and #20 in TABLE 
that can realize a gapped Z 2 spin liquid with up to next 
nearest neighbor mean-field amplitudes. 

Since the gapped spin liquid appears in p roximity to 
two magnetic orders in DMRG studies^ESl, it’s desir¬ 
able to understand the symmetry-breaking phases in the 
neighborhood of symmetric Z 2 spin liquids. In particular, 
a Z 2 spin liquid can enter either a noncollinear magn etic 
ordei*2212ni or a valence bond solid (VBS) phasd^^*^ via 
continuous quantum phase transitions, by condensing ei¬ 
ther bosonic spinons b or visons v. In this work, we study 
the magnetic and VBS phases in the neighborhood of the 
two gapped Z2 spin liquids, #1 and #20 in TABLE |Tj In 
particular we And that #20 can be driven into 120-degree 
noncollinear Neel order by a continuous quantum phase 
transition. 

In DMRG studies of triangular J 1 -J 2 model, strong ev¬ 
idence for nematic order i.e. brea king of 6-fold rotational 
symmetry has been observecP^. In this work we also ex¬ 
plore the possibility of gapped nematic spin liquidd^ in 
the neighborhood of symmetric spin liquids. In addition 
to the 2 symmetric candidates (#1 and #20), we iden¬ 
tify another promising nematic gapped state i.e. #6 in 
TABLE The mean-field ansatz of #6 and #20 nematic 
states are depicted in FIG. and 

This paper is organized as follows. In section |ll] we 
classify 20 distinct symmetric Z 2 spin liquids in the 
Abrikosov-fermion representation, as summarized in TA¬ 
BLE [T] We also unify the 8 previously obtained Z2 spin 
liquids in the Schwinger-boson representatioiJ^ with our 
Abrikosov-fermion states. In section |HI| we propose 2 
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promising candidates (^^1 and ^^20 in TABLE |T| among 
all 20 symmetric states for the J 1 -J 2 spin liquid phase, 
and study the magnetic orders and VBS patterns in the 
neighborhood of these 2 promising candidates. In section 
|IV| we introduce nematic order into symmetric states, to 
achieve nematic Zi spin liquids on triangular lattice. We 
find 3 promising gapped nematic Zi spin liquids ( 7 ^ 1 , 7^6 
and #20 in TABLEsee FIG. i an d§ that may be real¬ 
ized in triangular-lattice J 1 -J 2 model. Finally concluding 
remarks are given in section [Vj 

II. CLASSIFYING SYMMETRIC Z2 SPIN 
LIQUIDS ON TRIANGULAR LATTICE 


Therefore each physical symmetry operation u is fol¬ 
lowed by a SU{2) gauge rotation Guiv) when acting on 
fermionic spinons. It follows that distinct Z 2 spin liq¬ 
uids are classified by so-called projective symmetry group 
(PSG)I^ i.e. gauge-inequivalent symmetry operations 
{Gii(r)|it e symmetry group} with Ge(r) = ±1 G ^2 
for the identity element e of symmetry group. 

In our case of triangular lattice, the symmetry group 
(including spin rotations, time reversal T and space 
group generated by Ti 2 ,cr,G 6 ) imposes the algebraic 
conditions summarized on the left column of TABLE |ll] 
for symmetry operations {Gii(r) G SU{2)}. The gauge- 
inequivalent solutions are the following: 


We label each lattice site on a triangular lattice by 
its positional vector r = xai + ya 2 = ix,y), where 
ai = a(l,0) and 0,2 = a(—1,-\/3)/2 are chosen to be two 
Bravais vectors as shown in FIG.[2 The crystal symme¬ 
try group is generated by two translations Ti^ 2 , mirror 
reflection cr and site-centered 7 r /3 rotation Gg 

ix,y) ^ {x + l,y), {x,y) ^ {x,y+ 1), 

{x,y){y,x), {x,y) {x-y,x). (1) 

These crystal symmetries, together with global time re¬ 
versal T and SU{2) spin rotational symmetries impose 
constraints on a gapped Z 2 spin liquids. These symme¬ 
try conditions give rise to distinct Z 2 spin liquids in the 
Abrikosov-fermion representation. 


A. Classifying Z2 spin liquids in Abrikosov-fermion 
representation 


GT^{x,y) = 1 , GrAx^y) = r]f2, 

G^{x,y) = Gcs{x,y) = gce, 

GT{x,y)=gT- (5) 

with 

9l=VcT, gce = VCe, {gc^g^f = VcrC^, ( 6 ) 

5 t = VT, ga-gr = gTga-Va-T, gCegr = gTgCeVCeT- 

Notice that we’ve defined the time reversal action Gt{j^) 
on fermionic spinons f^.a in the following waj^^ 




GT(l')'r^$rO’*^, ‘hr = 


/r,T fr,i 

# _ ft 

Jr,I Jr,t 


In this notation, the mean-field amplitude (ri|r 2 ) trans¬ 
forms as 


(ri|r2) = -GT(ri)(ri|r2)G^(r2). (7) 


In Abrikosov-fermion representationP^^^, a spin 1 /2 on 
lattice site r is written in terms of two fermions /r.t/i- 



where a stands for Pauli matrices. First we construct 
mean-field Hamiltonian (“ansatz”) for fermionic spinons 

HMF = Er,,^^{n{rmi), (r|l> = (l|r)t. (3) 

where mean-field amplitudes (r|l) preserve all symme¬ 
tries. Then the many-spin wavefunction can be obtained 
by Gutzwiller projection which enforces the following 
single-occupancy constraint 

fr.tfr,1- + fl,lfr,i = 1 (4) 

on the mean-field groundstate \MF) of ferm ionic spinons. 
There is a SU{2) gauge redundancjPTlllISIl fj^jg repre¬ 
sentation: spin operators Sj. remains invariant under lo¬ 
cal SU{2) gauge rotation 'hr —>■ Wr^'r, Wr G SU{2). 


under anti-unitary time reversal T and 

{gri\gr 2 ) = Gg{gri){ri\r 2 )Gl{gr 2 ). (8) 

under any (global and spatial) unitary symmetry g. 

If r/T — ^ we have gx — I, and hence all mean-field am¬ 
plitudes must vanish identically due to time reversal sym¬ 
metry ([^. Therefore the only physical choice is rjT = —I 
and we can choose a gauge so that gx = ir^ = Gx{x, y), 
which only allows real hopping (r^) and real pairing (r^) 
amplitudes. With this we summarize all symmetric PSGs 
in TABLE They correspond to 20 distinct symmetric 
Z 2 spin liquids on triangular lattice. 

In addition to all symmetric states summarized in TA¬ 
BLE 1^ there is an extra pair of solutions to the PSG 
equations (|^ with 

ST = g^ = iT^, gcs = e'^^^ {ne = 1,2). 

as noticed in Ref. [5^ In particular the d-l- id paired state 
of Abrikosov fermions proposed in Ref. |T5] corresponds to 
the case of rig = 2 . 
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Symmetric Abrikosov-fermion states 

Nematic states: mean-field amplitudes 

Schwinger-boson states 

Label 

912 

Qcr 

gce 

onsite [0, 0] 

NN [1,1] 

NNN [2,1] 

NN [1,0] 

NN [1,1] 

NNN [2,1] 

NNN [1,-1] 

(P1,P2,P3) 

#1 

1 

r° 


^1.3 

^1,3 

^1,3 

^1.3 

^1,3 

^1.3 

fy,3 

(1.1,0) 

#2 

-1 

r° 

r° 

^1.3 

0 

0 

0 

0 

0 

0 

(0,1,0) 

#3 

1 

r° 


0 

0 

0 

0 

0 

0 

0 


#4 

-1 

r° 


0 

0 

^1,3 

^1.3 

0 

^1.3 

^1,3 


#5 

1 

r° 

ify 








(1.1,1) 

#6 

-1 

r° 

ify 


0 



0 



(0,1,1) 

#7 

1 

ir^ 

r° 

0 

0 

0 

^1.3 

0 

^1,3 

0 


#8 

-1 

ir^ 

rO 

0 

0 

0 

0 

0 

0 

0 


#9 

1 

ir^ 


0 

0 

0 

0 

0 

0 

0 


#10 

-1 

ir^ 


0 

^1,3 

0 

^1.3 

^1.3 

^1.3 

0 


#11 

1 

ir^ 


0 

0 

0 


0 


0 


#12 

-1 

ir^ 

ify 

0 


0 




0 


#13 

1 

1# 

r° 




^1.3 

r® 

^1,3 


(1,0,0) 

#14 

-1 

i# 

tO 


0 

0 

0 

0 

0 

0 

(0,0,0) 

#15 

1 

1# 


0 

0 

0 


0 


0 


#16 

-1 

1# 


0 

0 



0 




#17 

1 

i# 


0 

0 

0 

0 

0 

0 

0 


#18 

-1 

i# 


0 



^1.3 


^1.3 



#19 

1 

1# 






r® 



(1.0,1) 

#20 

-1 

i# 

ify 



0 




0 

(0,0,1) 


Table I: Summary of 20 distinct symmetric Z 2 spin liquids on triangular lattice. Among them, 8 states have their counterparts 
in the Schwinger-boson representatiorP^. We use \x,y\ = {0,0|a;,j/) to label representative mean-field amplitudes. Only 4 
Abrikosov fermion states support a symmetric Z 2 spin liquid with up to next nearest neighbor (NNN) mean-field amplitudes, 
as highlighted by red color. Each symmetric Abrikosov-fermion state have a unique nematic descendant which breaks Ce 
rotation but preserves mirror reflections cr and R = rr(C6)^. In addition to the 4 symmetric ansatz, another 4 states support a 
nematic Z 2 spin liquid with up to NNN amplitudes, as highlighted by blue color. In the mean-field amplitudes, represents 
real pairing while stands for real hopping. Up to NNN mean-field amplitudes, only 3 ansatz support a gapped nematic Z 2 
spin liquids: they are ^1, #6 and ^^20 states. 


Once the symmetry transformations of fermionic 
spinons are determined, the associated mean-field am¬ 
plitudes [x,y\ = (0,0|a:, j/) must obey certain symme¬ 
try constraints. For example, on-site chemical potential 
/t = [0, 0] must satisfy 

5tA5t = -A, (9) 

A = S-rAffl = SCeA^Ce- 

Mean-field amplitude a = [1,1] between nearest neigh¬ 
bors (NNs) must have 

9 Tag\. = -a, (10) 

a = r]i2gcTagi, 

= ??l25Ce«5c6 ■ 

For next nearest neighbor (NNN) mean-field amplitude 
jS = [2,1] we have 

5 t/ 35 t = ( 11 ) 

m2gCeP9Ce = 

/ 3 I = gi29cj9ce- 


Among the 20 different Abrikosov fermion states sum¬ 
marized in TABLE |T] only 4 states can realize a sym¬ 
metric Z 2 spin liquid with up to NNN mean-field ampli¬ 
tudes. They are fyl, fy6, fyl8 and fy20 as highlighted by 
red color in TABLE |T] Later we’ll show that fyl8 state 
doesn’t support a gapped spinon spectrum with up to 
NNN amplitudes. 


B. Duality between Abrikosov-fermion and 
Schwinger-boson representations 

In Schwinger-boson represent atioiP^, a spin-1/2 par¬ 
ticle on lattice site r is decomposed into two species of 
bosonic spinons {fy^aja =t / i}: 

Sr = 2 ^ (12) 

a,/3=t/4 

where a are Pauli matrices. Similar to Abrikosov- 
fermion representation, the spin wavefunction is obtained 
by single-occupancy Gutzwiller projection on mean-field 
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Figure 1: (color online) Crystal symmetry of triangular 
lattice, and nearest-enighbor (NN) hopping amplitudes of 
fermionic spinons in the tt-Aux state on the triangle lattice. 
All NN hoppings are real with the same magnitude |q| > 0, 
where thick (thin) bonds have negative (positive) hopping 
signs. Each magnetic unit cell (denoted by the green oval) 
is chosen to contain 4 lattice sites labeled by red/blue dia¬ 
monds/dots. 


states of bosonic spinons. Schwinger-boson representa¬ 
tion has a U{1) gauge redundancy e'‘^‘'6r,cr and 

hence any symmetry operation u is followed by a C/(l) 
phase rotation of Schwinger bosons by angle (j)u{x,y). 
Ref. I3H classified Z 2 spin liquids in Schwinger-boson 
representation and obtain 8 different states on triangu¬ 
lar lattice, with the following gauge-inequivalent PSGs 
fe = 0 ,l): 

{x, y) = 0, (j)T 2 {x, y) = xpiTT, (13) 

(l>(T{x,y) = {xypi +P 2 / 2 )tt, 

<Pce{x,y)= [piy{x+^)+p3/2]TT. 

and one can always choose a gauge so that 4>T{x,y) = 0 
for time reversal T. 

What is the relation between these 8 Schwinger-boson 
states and the 20 Abrikosov-fermion states obtained in 
this work? There are two ways to achieve this connec¬ 
tion. The first approach is to directly compare the crystal 
symmetry (mirror reflection cr^R and inversion (Ce)^) 
quantum numbers of their projected spin wavefu nction 
on different finite lattice^^ and on infinite cylinder^2llll! 
For symmetry quantum numbers associated with time re¬ 


Algebraic Identities 

bosonic ba 

fermionic fa 

vison V = b X f 

T-It-It^Ti = e 


rii2 

-1 

cr~^TicrTp^ — e 

1 

1 

1 

cr-^T 2 crT{-^ = e 

1 

1 

1 

Ce ^T\CqT 2 = e 

1 

1 

1 

= e 

1 

1 

1 

= e 

(_1)P2 


1 

= (G6(t)= = e 


V^rCe 

1 

(Ge)® = e 

(_1)P3 

VCe 

-1 

Tf iT"iriT = e 

1 

1 

1 

T-^T~^T 2 T = e 

1 

1 

1 

o-"1T“VT = e 

(_1)P2 

PerT 

1 

R-^T-^RT = e 

^_]^^P2+P3 

VCeTPcrT 

1 

= e 

-1 

-1 

1 


Table II: The correspondence between bosonic spinorP^, 
fermionic spinon and vison PSGs on the triangular lattice. 
Here bosonic spinon {pa) PSGs are labeled by three integerP^ 
Pi = 0,1 (i = 1,2,3), while fermionic spinon {fa) PSGs are 
labeled by six integers (1712, Per, PcrCg, ? 7 C 6 , ho-T, PCst) where 
p = ±1. Ghoosing a proper gauge we can always fix 
C 6 ~^T 2 CeTp^Tf~^ = Ce~^TiCeT 2 = 1 for both spinons and 
visons. 


versal symmetry (3 row at the bottom of TABLE [ll|, one 
can compare the degeneracy in the entanglement spec¬ 
tra of their projected wavefunctions on different infinite 
cylinders with cr and R symmetrie^^. Repeating the 
same analysis as in Ref. | 3 S] we obtain the unification of 
Schwinger-boson PSGs and Abrikosov-fermion PSGs as 
shown in TABLE HI 

The 2nd approach is to understand the vison PSGs 
in the Schwinger-boson state. As argued in Ref. |40] 
any Z 2 spin liquids constructed by projecting a mean- 
field Schwinger-boson state doesn’t support gapless edge 
stated^. On the other hand, nontrivial vison PSGs, 
e.g. (T^ — —1 or crTcr~^T~^ — —1 acting on visons u, 
will lead to symmetry protected gapless edge states on 
an open edge preserving reflection cr. These strong con¬ 
straints can determine the vison PSGs in any Schwinger- 
boson state, as shown in the right column of TABLE [lT| 
Meanwhile, in a Z 2 spin liquid the fermionic spinon 
can be viewed as a bound stat^i^ of a bosonic spinon 
ha and vison v due to fusion rule fa = ba v. Gon- 
sequently the fermionic spinon (/) PSGs is a product 
of bosonic spinon {b) P SGs and vison {v) PSGs, mul¬ 
tiplied by a twist factoii^SIllHini, This twist factor is - 
1 in the case of (Ge)® = cr'^ = R^ = #SH40l 
crTcr-iT"^ = RTR-^T~^ = and +1 in all other 

cases. 

As a result in TABLE [IT] a Schwinger-boson state 
{Pi,P 2 tP 3 ) belong to the same symmetric Z 2 spin liquid 
phase as an Abrikosov-fermion state if their PSGs have 
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the following correspondence: 

^i2 = (-ir, (14) 

= V^T = VCellcrCs = 

VCe = VCeT = ■ 

All 8 Schwinger-boson states have their counterparts in 
Abrikosov-fermion representation, as summarized in TA¬ 
BLE lU Since the continuous phase transitions between 
Z 2 spin liquids and noncollinear magnetic order s can be 
simply achieved by condensing bosonic spinond^siMI 
the Schwinger-boson representation, this unification al¬ 
lows us to understand the neighboring magnetic orders 
of some Abrikosov-fermion states. 

C. Nematic descendants 

DMRG studied^E^ suggest strong evidence for nematic 
order in the gapped spin liquid phases, i.e. spontaneously 
broken Cq rotational symmetry coexisting with Z 2 topo¬ 
logical order. This motivates us to study nematic (or 
anisotropic) Z 2 spin liquids on a triangular lattice. In 
spite of broken Cq symmetry, these nem atic Z 2 spin 
liquids still preserve two mirror reflectioni j^ l ^° ( cr and 
R = as well as their combination: site-centered 

inversion symmetry I = (Cq)^. Ref. classified ne¬ 
matic Z 2 spin liquids on an anisotropic triangular lattice 
(with (T and R symmetry), leading to at least 63 distinct 
Abrikosov-fermion states. Instead of examining all these 
states, here we focus on those nematic states in proximity 
to a Ce-symmetric Z 2 spin liquids. In other words they 
can be achieved by adding anisotropic perturbations to 
a Cg-symmetric Z 2 spin liquid, i.e. one of the 20 states 
in TABLE n 

It turns out each of the 20 Abrikosov-fermion states 
summarized in TABLE [I] has a unique nematic descen¬ 
dant with the following symmetry transformations 

GtA^^v) = GTAx,y) = ■ni2, 

G^{x,y) = Gi{x,y) = Vi2^{gcJ^, 

Gr{x, y) = gn, Gt{x, y) = gr = ir^. 

where gu = ga-igce)^- These symmetry properties im¬ 
pose the following constraints on a generic mean-field 
amplitude [x,y] = {0,0\x,y): 

gT[x,y]gT = - [x,y], 

g^[x,x]gl = <2 [a:, a:], 

9R[x,-x]g]^ = gf2{x,-x]. (15) 

To be specific, there are two independent amplitudes be¬ 
tween 1st nearest neighbors (NNs) i.e. [1,0] and [1,1] in 
nematic states. They satisfy 

5T[l,0]5i. = -[l,0], 

9UG0]9c^ = m2[G0Y. ( 16 ) 


and 

9t[1,1]9t = -[1)1]. 
ffCe[l>l]ffC6 = '7l2[l.l]^ 

5<T[l,l]gi = ??12[1,1]- (17) 

Similarly, two independent 2nd NN amplitudes [2,1] and 
[1,-1] satisfy 

gT[2,l]ff; = -[2,l], 

9hi^A]9c^ = 91212, IV. (18) 

and 

9t[1,-1]9t = -[l.-lj. 

9c^^,-^]9c! = 9i2[i,-lV, 

5.[l,-l]fft = [l,-l]t. (19) 

These conditions give rise to symmetry-allowed mean- 
field amplitudes for the 20 nematic Z 2 spin liquids, as 
summarized in TABLE Due to spontaneous breaking 
of Ce rotational symmetry, the symmetry conditions for 
nematic Z 2 spin liquids are less stringent than for the 
isotropic (Cg symmetric) ones, leading to nonzero NN 
and NNN amplitudes that are not allowed in isotropic 
states. As a consequence, 4 more states support a ne¬ 
matic Z 2 spin liquid with up to NNN mean-field ampli¬ 
tudes, in addition to the 4 isotropic states (red-colored 
in TABLE |T| in isotropic case. They are #4, ^7, #10 
and #13, as highlighted by blue color in TABLE |T] 

Previous variational Monte Carlo studjEl suggested a 
gapless d-wave paired state of Abrikosov fermions, which 
is energetically competitive for a finite range of J 2 /J 1 . 
In this state the fermionic spinons share the same PSG 
with state #13 in TABLE |T] 

III. GAPPED SYMMETRIC Z 2 SPIN LIQUIDS 
AND THEIR NEIGHBORING PHASES 

In DMRG studies of spin-1/2 J 1 -J 2 Heisenberg model, 
the gap ped s pin liquid phase appear only in the parame¬ 
ter rang^^ESlof J 2 /J 1 < 0.17. Therefore for candidate Z 2 
spin liquids in Abrikosov-fermion representation, we fo¬ 
cus on Z 2 spin liquid states that can be realized by up to 
next nearest neighbor (NNN) mean-field amplitudes. In 
particular we’ll start from mean-field ansatz with nonzero 
nearest neighbor (NN) amplitudes. We can always choose 
a proper gauge such that NN mean-field amplitudes are 
all real hoppings (r^). 

There are only two possible NN hopping ansatz that 
preserve all symmetries (time reversal, spin rotation and 
crystal symmetries): 

(I) uniform RVB state: all NN real hoppings share the 
same amplitude a. 

(II) TT-flux (algebraic spin liquid) state: tt flux is in¬ 
serted in half of the triangles (see FIG. [^. 
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They correspond to two distinct symmetric U{1) spin 
liquids on triangular lattice, which host gapless fermionic 
spinons coupled to t/(l) gauge field. In the following we’ll 
study candidate Z 2 spin liquids by adding perturbations 
to these two parent U{1) spin liquids. 


A. Gapped Z2 spin liquids near uniform RVB state 


Among the 20 states in TABLE 4 of them support 
uniform NN real hoppings, corresponding to 4 symmet¬ 
ric Z 2 spin liquids in the neighborhood of uniform RVB 
state. They are #1, #5, #13 and #19 states. Three of 
these 4 states (#5, #13 and #19) only allow real hop¬ 
pings up to NNN amplitudes, therefore remaining to be 
gapless U{1) spin liquids with spinon fermi surfaces. 

On the other hand, #1 state allow uniform hoppings 
and pairings in NN and NNN amplitudes. It’s straight¬ 
forward to verify that uniform pairing terms can open up 
a gap on the spinon fermi surface of NN hoppings, giving 
rise to a gapped Z 2 spin liquid which preserves all sym¬ 
metries. Therefore #1 state, which is dual to Schwinger- 
boson (1,1,0) state, is the only gapped Z 2 spin liquid 
phase in the neighborhood of uniform RVB state. 


B. Gapped Z2 spin liquids near tt-Aux algebraic 
spin liquid 


Unlike the uniform RVB state with a spinon fer mi sur- 
face, the tt-Aux statd^ is an algebraic spin liquicP^*^, 
described by gapless Dirac fermions coupled to emergent 
t/(l) gauge field. Similar to the case of uniform RVB 
state, there are 4 distinct Z 2 spin liquids in the neigh¬ 
borhood of TT-flux state: i.e. #10, #12, #18 and #20 
states in TABLE HI 

Choosing the 4-site magnetic unit cell as illustrated in 
FIG.[^ the Bloch Hamiltonian in the basis of '4>{ri,r2),<y — 

(/*(2x,2y-t-l),(T ) /'(23:-t-l,2y-t-l),cn f{2x^2y),aj /*(2x-t-l,2y),(T) 

writes 


HIif = a Ek.a K = 






gifc2 _ g“ifci 

1 + e 

i k\ 

0 

gi(/ci-|-fc2) _ 2^ 

1 

1 + e“ 

- i k2 

g-i(fci-|-fc2) _ ^ 

0 

_l_g-ifci 

\^e-ife2 _ 

gifci 


-1-gifci 

0 


The dispersion crosses zero energy (fermi level) at mo¬ 
mentum fci = k 2 = TT. Expanding the above Bloch 
Hamiltonian around Dirac point (tt, tt) leads to the 
following linearized Dirac Hamiltonian in the basis of 


^q+( TTjTt) ,<t 

hlc= [ - qifJ-zVy + q2^J■yV0 + (?! + 92)Mxt'j/] CQ + ^dq^) 

= '/&{qxix + qyly) + 0(|qp), 

_ (yi:-2.n^)vy-yyVo _ y:,:Vy+y.yVo 

Ve ’ ly - U2 ’ 
q = k - (7r,7r) = ( 91 , 92 ) = {2qy„y/lqy - q^,). 

where v and /2 are Pauli matrices for the (diamond,dot) 
and (red,blue) sublattice indices respectively (see FIG. 
[^, while tf are Pauli matrices associated with spin index. 
We’ve set lattice constant as unity. 

In the Nambu basis of 4>q = — i#Lq^o-crj/)^, the 

low-energy Dirac Hamiltonian writes 

( 20 ) 

<iq = [ - ql^J^zl^y + 92M!/^0 + (9l + 92)Mxt'y] CToT^ 

= VQ[qxlx + qylyWoTz- 


where r are Pauli matrices for Nambu index. 

A gapped Z 2 spin liquid can be obtained by adding 
pairing terms, which open up a gap in the Dirac spec¬ 
trum. These pairing “mass terms” for Dirac fermions, 
however, may break certain symmetries. In particular, 
Dirac fermion 4'q = cTyTy'^'Lq in (201 transform under 
symmetries in the following way: 


4'q ^ icry4'-q, 

4'q —)■ exp(i|n • cr)4'q, 

4'q ^ MO^'yCToTj^'q, 

^( 91192 ) ^<tO’o(9<t)'I^( 132 , 91)1 

^( 91 , 92 ) 

'^( 91 , 92 ) ^ ^Ce^^oigCs)^(-q2,qi+q2)- 


where we define 9 ^. = and gcg = 

giTTT in a different gauge. We have 

{9Ce,'rz} = { 9 ct,Tz} = 0 and 


U^ = Ur 





0 

0 

A 



0 

0 


0 

-1 

0 

0 

TJ^ — 

1 

0 

0 

0 

0 

0 

1 

0 

5 ^ Cq 

0 

0 

1 

0 

Vi 

0 

0 

V 



-1 

0 

0/ 




Notice that there is a U{1) gauge redundancy for alge¬ 
braic spin liquid (20), i.e. the Dirac Hamiltonian is in¬ 
variant under any Tz rotations. 

Glearly the symmetric spin-singlet pairing mass terms 
which anticommute with Dirac matrices "fx.yO'oTz are sim- 

ply 


— I^O^O^oip^^ssc ‘ '^x “1“ '^^ssc ‘ (^1-) 
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Time-reversal-invariant real pairing mass pre¬ 

serves (T and Cg symmetries, if and only if [gceT'''x] = 
[gcr,Tx] = 0. Therefore among the 4 neighboring Z 2 
spin liquids near tt-Aux state, only ^20 state with g^- = 
gCa = It’z allows a symmetric pairing mass for the Dirac 
fermion. 

As a result, ^20 state which is dual to Schwinger-boson 
( 0 , 0 , 1 ) state is the only gapped Z 2 spin liquid in the 
neighborhood of the 7r-flux algebraic spin liquid. 


C. Confined symmetry-breaking phases in 
proximity to gapped Z2 spin liquids 

A gapped Z 2 spin liquid can be driven into “confined” 
symmetry-breaking phases with no fractional excitations, 
by condensing their fractional excitations. For example, 
condensing spin-singlet vison excitations typically gives 
rise to a valence bond solid (VBS) phase, which breaks 
crystal symmetries but preserves global spin rotation and 
time reversal symmetries. On the other hand, condensing 
bosonic spinons (carrying spin- 1/2 each) typically leads 
to a noncolinear magnetic order which breaks time rever¬ 
sal as well as spin rotational symmetries. Therefore con¬ 
fined symmetry-breaking states in proximity to a Z 2 spin 
liquid reflect the symmetry quantum numbers of spinons 
and visons therein, and is an important characterization 
of the symmetry fractionalization in the Z 2 spin liquid. 
In the following we analyze the proximate ordered phases 
of the two promising Z 2 spin liquids, 7/1 and /^20 as dis¬ 
cussed above. 

First we discuss 7/1 state in the neighborhood of uni¬ 
form RVB state. As shown in section |IIB[ it describes 
the same Z 2 spin liquid phase as the ( 1 , 1 , 0 ) state in 
Schwinger-boson represent atioiP^. In the Schwinger- 
boson representation, the transition from Z 2 spin liquid 
to a magnetic ordered phase is simply described by con¬ 
densation of Schwinger-bosons. As shown in Ref. [Ml the 
consequent ordered phase in proximity to ( 1 , 1 , 0 ) state is 
a 4-sublattice noncolinear magnetic order: 

S{r) = ni(-l)" -b n 2 (-l)*' + n 3 (-l)"+^ 

It features a 2 x 2 magnetic unit cell and is the tetrahedral 
state classified in Ref. US] 

On the other hand, the low-energy dynamics of spin¬ 
less visons is described by a fully frustr ated quantum 
Ising model on the dual honeycomb lattic^^Hm^ which is 
constrained by the vison PSGs in TABLE |ll] As shown 
in Ref. l32l44fH5l the VBS phase achieved by condensing 
visons in the Z 2 spin liquid is featured by an enlarged 
■\/l2 X y /12 unit cell. 

Now let’s turn to ^20 state in the neighborhood of 
TT-flux algebraic spin liquid. It belongs to the same 
Z 2 spin liqui d pha se as (001) state in Schwinger-boson 
representatiorpSlMI Since all Z 2 spin liquids constructed 
from Schwinger-boson mean-field ansatz share the same 


vison PSG^^, the VBS phase in proximity to ^20 state 
generically has the same VT2 x symmetry-breaking 
pattern as the /^1 state. Meanwhile upon condensation 
of bosonic spinons. Ref. I29I34I showed that the so-called 
120-degree Neel ordered phase with a ■\/3 x magnetic 
unit cell will be developed from the ( 001 ) state. 

Since ^20 state can be obtained by adding a singlet 
pairing mass term to Dirac spinons in tt-Aux state, the 
Abrikosov-fermion representation provides an alterna¬ 
tive way to understand the confined symmetry-breaking 
phases in its neighborhood. To be specific, the VBS 
phase corresponds to the following mass term added to 
Dirac spinons in (20 1 : 


Myss = B ■ V = ^ Bi{g.x -\- Hz — gyVy) + (22) 


B2{Vx: + gxVz - g-zVz) + + fJ-z^x) T, 


where we’ve only written down the Pauli matrices with a 
nonzero subscript. The VBS mass matrices V are invari¬ 
ant under time reversal and spin rotations, but transform 
nontrivially under space group symmetries: 

ri(Pl,V2,V3)Tfl = (Vu-V2,-V2), 

T2{Vi,V2, V 3 )T-^ = (-Vl, -V2, V2), 
(t{VuV2,V3)o--^ = (-V3,V2,-Vi), 
Ce{Vi,V2,V3)Ce ^ = {—V2,V3,—Vi). 

corresponding to crystal-symmetry-breaking VBS orders. 
It’s straightforward to verify that 3 VBS masses Vi, 2,3 
anticomute with each other, and they all anticommute 
with the singlet pairing mass Tx,y corresponding to Z 2 
spin liquid phase. As a result the low-energy dynamics 
of 5-component mass-term vector 


V = VBSi,VBS 2 , VBS 3 ) 


is described by 0(5) non-linear sigma model (NLsM) 
with a Wess-Zumino-Witten (WZW) terirP^HUl^ Such a 
WZW term can be derived by integrating out fermionic 
spinons in the Dirac Hamiltonian (|20t. It implies a con¬ 
tinuous quantum phase transitiorl^® 2 l between Z 2 spin 
liquid #20 (A^sc ^ 0, VBS = 0) and the VBS phase 
(A,,, = 0, VBS^Q). 

Similarly the quantum spin Hall (QSH) mass terms to 
Dirac spinons 


Mgs// = n • m = 


iHxJ^y + g-zVy - gy){n ■ a)Tz 




(23) 


also anticommute with the singlet pairing masses. This 
indicates the dynamics of the following 5-component 
mass vector 


v' = {^Assc, QSHx, QSHy, QSHz) 

is also captured by 0(5) NLsM with a WZW term. Once 
the Dirac spinons acquires a finite gap from the QSH 
mass term, the monopole of U{\) gauge fields coupled 



















Figure 2: (color online) Mean-field hopping amplitudes (up to 
NNN) of nematic Z 2 spin liquid #20 in TABLE There are 
two types of NN hoppings denoted by black and red colors: for 
black ones, thick and thin lines have the same amplitude |ai| 
but opposite signs; for red ones, solid and dashed lines have 
the same amplitude |q 2 | but opposite signs. NNN amplitudes 
are labeled by blue color: solid and dashed lines have the 
same amplitude |/3i| but opposite sings. Onsite real paring is 
also allowed by symmetry. 


Figure 3: (color online) Mean-field hopping amplitudes (up 
to NNN) of nematic Z 2 spin liquid #6 in TABLE Nonzero 
NN hoppings are denoted by black thick lines: for black ones, 
solid and dashed lines have the same amplitude |ai| but oppo¬ 
site signs. The thin lines connecting NNs have zero hopping 
amplitude. There are two types of NNN amplitudes: for blue 
(or red) colors, solid and dashed lines have the same ampli¬ 
tude |/3i| (or I/I 2 I) but opposite sings. Onsite real paring is 
also allowed by symmetry. 


to fermionic spin ons will carry a spin quantum number 
due to QSH effectl^^l^^*^. Although QSH mass Mgsjy in 
(231 appear to preserve C/(l) spin rotation along n-axis, 
the proliferation of monopole events in U{1) gauge fluc¬ 
tuations will destroy the spin conservation along h-axis, 
leadin g to a noncolinear magnetic order with 3 Goldstone 
modeJ^^*^. Therefore the WZW term for the quituplet 
of singlet pairing and QSH masses implies a continuous 
quantum phase transition between 120-degree Neel order 
and Z 2 spin liqu id #2 0, in accordance with Schwinger- 
boson (001) stat^2llS^ 


In the phase diagram of triangular J 1 -J 2 Heisenberg 
model obtained in DMRG studies, the gapped spin liq¬ 
uid phase appear in proximity to 120-d egre e Neel order 
and another colinear “stripy” Neel order*21i^. This is an¬ 
other evidence in favor of #20 state {i. e. Schwinger-boson 
(0,0,1) state) as a promising candidate for the spin liquid 
found in J 1 -J 2 model. 


IV. NEMATIC ^2 SPIN LIQUIDS ON 
TRIANGULAR LATTICE 

In DMRG studies of spin-1/2 J 1 -J 2 Heisenberg model 
on triangular lattice, both the even and odd sectors of 
the gapped spin liqu id p hase show spatially anisotropic 
spin-spin correlations^iSl^ different from the Z 2 spin liq¬ 
uid in kagome Heisenberg modeP. Meanwhile the spin 
liquid shows strong response to a Gg-breaking perturba¬ 


tion to Heisenberg couplingP. This motivates us to seek 
for nematic Z 2 spin liquid candidates for J 1 -J 2 Heisen¬ 
berg model. Again we focus on those nematic Z 2 spin 
liquids that can be realized by up to NNN mean-field 
amplitudes. 

Nematic Z 2 spin liquids on triangular lattice has been 
classified in Ref. [33l Here we restrict ourselves to those 
nematic (anisotropic) Z 2 spin liquids that can be ob¬ 
tained by perturbing fully-symmetric (isotropic) Z 2 spin 
liquids. As discussed in section |H G| for each symmet¬ 
ric Z 2 spin liquids on triangular lattice, there is one and 
only one nematic “descendant” which preserves two mir¬ 
ror reflections cr and R but breaks Cq rotation (see FIG. 
[^. Their symmetry-allowed mean-field amplitudes are 
summarized in TABLE HI 

First we examine those nematic Z 2 spin liquids in the 
neighborhood of uniform RVB state, #1, #5, #13 and 
#19. In addition to #1 state which already realizes a 
gapped Z 2 spin liquid without breaking Cq symmetry, in 
the presence of nematic order #13 state can also real¬ 
izes a Z 2 spin liquid with up to NNN amplitudes. To be 
precise, Cq symmetry breaking gives rise to nonzero pair¬ 
ing amplitudes (t^) in #13 state, as shown in TABLE 
However, these nematic pairing terms generally cannot 
open a full gap on the spinon fermi surface: typically 
they leave 4 point nodes on the fermi surface similar to 
the d-wave superconductors. Therefore even in the pres¬ 
ence of nematic order, #1 state is still the only gapped 
Z 2 spin liquid near uniform RVB state. 

Next we look into nematic Z 2 spin liquids in the 
neighborhood of tt-Aux state, #10, #12, #18 and 
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^20. Among them, only ansatz ^20 supports a fully- 
symmetric gapped Zi spin liquid. After breaking 
symmetry with nematic order, it is possible to realize 
a Z 2 spin liquid with up to NNN amplitudes in both 
7^10 and 7^18 (but not #12) state. However, the pairing 
terms cannot open up a gap in the spinon Bogoliubov 
spectrum: point nodes still exist in #10 and #18 state. 
As a result, #20 is again the only gapped nematic Z-i 
spin liquid near tt-Aux algebraic spin liquid. In addition 
to onsite real pairing, it also allows NN/NNN hoppings 
as demonstrated in FIG. [21 

Meanwhile, there are 3 other nematic Zi spin liquids 
that can be realized with up to NNN amplitudes: #4, 
#6 and #7. Among them, #4 and #7 have gapless 
excitations in their spinon Bogoliubov spectra, while #6 
supports a gapped nematic Zi spin liquid. The mean- 
field ansatz #6 allows onsite real hoppings, together 
with NN/NNN real hoppings as demonstrated in FIG.|^ 

To summarize, there are 3 gapped nematic Zi spin 
liquids that can be realized by up to NNN mean-field 
amplitudes among the 20 states in TABLE [^ they are 
#1 state dual to Schwinger-boson (110) state, #6 state, 
and #20 state dual to Schwinger-boson (001) state. They 
are the most promising candidate for possible nematic Zi 
spin liquid phase in triangular-lattice J 1 -J 2 Heisenberg 
model. 


V. SUMMARY 

To summarize, in this work we studied symmetry 
Z 2 spin liquids on triangular lattice and their nematic 
descendants, in connection to the numeric discovery 
of a gapped spin liquid phase in triangular lattice 
J 1 -J 2 model. We classified symmetric Z 2 spin liquids 
in the Abrikosov-fermion representation, found 20 
distinct states and identified 2 states (#1 and #20 in 
TABLE |T]) as the promising candidates which support 


a gapped Z 2 spin liquid with up to NNN mean-field 
amplitudes. By establishing the correspondence between 
symmetric Abrikosov-fermion and Schwinger-boson 
states on triangular lattice, we are able to understand 
the noncollinear magnetic order phases and VBS phases 
in the neighborhood of these 2 gapped spin liquids. In 
particular state #20 in TABLE dual to the (001) 
state in Schwinger-boson representation, is continuously 
connected to the 120-degree Neel order, which is found 
to lie in proxi mity to gapped spin liquid in the DMRG 
phase diagrairplifli. Motivated by numerical evid ence 
of possible nematic order in this spin liquid phas^^E^ 
we also studied nematic Z 2 spin liquids by perturbing 
the symmetric ones. We identified 3 promising nematic 
states, #1, #6 and #20 in TABLE [^ which can realize 
a gapped Z 2 spin liquid with up to NNN amplitudes. 
It’ll be interesting for future variational Monte Garlo 
studies to see whether these nematic Z 2 spin liquids are 

energetically more favorable than fully-symmetric states 
or not. 

Note added: Upon completion of this paper, we be¬ 
came aware of an independent worlP^ which also classi¬ 
fied symmetric Z 2 spin liquids in the Abrikosov-fermion 
representation on triangular lattice. 
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